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The directed polymer in a 1 + 3 dimensional random medium is known to present a disorder- 
induced phase transition. For a polymer of length L, the high temperature phase is characterized 
by a diffusive behavior for the end-point displacement R 2 ~ L and by free-energy fluctuations of 
order AF(L) ~ 0(1). The low-temperature phase is characterized by an anomalous wandering 
exponent R 2 / L ~ L u and by free-energy fluctuations of order AF(L) ~ where u ~ 0.18. In this 
paper, we first study the scaling behavior of various properties to localize the critical temperature T c . 
Our results concerning R 2 /L and AF(L) point towards 0.76 < T c < T 2 = 0.79, so our conclusion 
is that T c is equal or very close to the upper bound T2 derived by Derrida and coworkers (T2 
corresponds to the temperature above which the ratio Zj j /(Zi,) 2 remains finite as L — » 00). We 
then present histograms for the free-energy, energy and entropy over disorder samples. For T> r c , 
the free-energy distribution is found to be Gaussian. For T <C T c , the free-energy distribution 
coincides with the ground state energy distribution, in agreement with the zero-temperature fixed 
point picture. Moreover the entropy fluctuations are of order AS ~ L 1 ^ 2 and follow a Gaussian 
distribution, in agreement with the droplet predictions, where the free-energy term AF ~ is a 
near cancellation of energy and entropy contributions of order L 1 ^ 2 . 

I. INTRODUCTION 

The model of a directed polymer in a random medium plays the role of a spin glass toy model in the field of 
disordered systems Q, 0, Sfl @- At low temperature, there exists a disorder dominated phase, where the order 
parameter is an 'overlap' @jEIl4 0- ■"■ n finite dimensions, a scaling droplet theory was proposed m direct 

correspondence with the droplet theory of spin-glasses whereas in the mean-field version of the model on the 
Cayley, a freezing transition very similar to the one occurring in the Random Energy Model was found |2j ■ The phase 
diagram as a function of space dimension d is the following p|. In dimension d < 2, there is no free phase, i.e. any 
initial disorder drives the polymer into the strong disorder phase, whereas for d > 2, there exists a phase transition 
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between the low temperature disorder dominated phase and a free phase at high temperature 0, [n| . This phase 
transition has been studied exactly on a Cayley tree B and on hierarchical lattice In finite dimensions, bounds on 
the critical temperature T c have been derived [HEaE! : T o(d) <T C < T 2 (d). The upper bound T 2 (d) corresponds 
to the temperature above which the ratio Z£/(Zl) 2 remains finite as L — > 00. The lower bound To corresponds to 
the temperature below which the annealed entropy becomes negative. On the Cayley tree, the critical temperature T c 
coincides with To 0. In finite dimensions, there is a debate to know whether T c coincides or not with T 2 . Arguments 
in favor of the equality T c — T 2 in finite dimensions d where ui(d) > can be found in [l5|. whereas a new upper 
*^ ■ bound T* based on the entropy of random walks was recently proposed in (l6|. As explained in detail in [l7|. the 
debate actually concerns the form of the negative tail of the free-energy distribution in the high temperature phase. 
If this tail decays only exponentially, as on the Cayley tree, then T c < T 2 is possible; if it decays more rapidly than a 
simple exponential one gets T c = T 2 . 

In this paper, we study numerically the scaling properties of the directed polymer in dimension 1 + d with d = 3. As 
for previous numerical studies of the freezing transition in d = 3, we are o nly aware of ref. |l3ll2l| . which draw different 
conclusions for the specific heat exponent a, namely a ~ —0.1 ± 0.1 in [y| and a = 2 — v ~ —2 ± 0.7 in .21]. With 
our numerical data, we are not able to give a reasonable quantitative estimate of this critical exponent (see below). 
However, our results concerning the scaling of the transverse displacement R 2 (L) and of the free-energy fluctuations 
AF(L) point towards a critical temperature T c which is equal or very close to the upper bound T 2 discussed above. We 
also present histograms of free-energy, energy and entropy as a function of temperature, in order to better understand 
the nature of the high and low temperature phases. 

The paper is organized as follows. In Sect [HI we define the model and give some numerical details. In Sect. IIIII 
we present our results concerning the end-point fluctuations. In Sect. IIVI we discuss the scaling behavior and the 
probability distribution of the free-energy, the energy and the entropy as T varies. In Sect. [V] we present the finite-size 
scaling analysis of spatial and thermodynamic observables. We summarize our conclusions in Sect. IVII 
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II. MODEL AND NUMERICAL DETAILS 

In this paper, we present numerical results for the random bond version of the model defined by the recursion 
relation on a cubic lattice in d = 3 

2d 

Zl+i{?) = J2 e~^ {P+S ^ rl Z L (f+ Sj) (1) 
The bond energies ei(r + ej,r) are random independent variables drawn from the Gaussian distribution 



The numerical values of the bounds discussed in the Introduction are then given by |13| 

T a {d = 3) = 0.528.. <T C < T 2 (d = 3) = 0.790... (3) 

In the previous numerical study of the same model |l3j , a finite size scaling analysis with respect to the transverse 
direction gave a critical temperature of order T c ~ 0.6, whereas our data presented below point towards T c ~ Ti = 0.79. 
In the high temperature phase (T > T c ), the free-energy density 

/(T)E-Tlimi^ (4) 

is known to coincide with the annealed value 

/CO = f ann (T) = -Tln(2d) - for T >T C (5) 

As a consequence, the energy e(T) and the specific heat c(T) also coincide with their annealed values e ann (T) = — 1/T 
and c ann (T) = 1/T 2 , so that the regular parts of thermodynamic quanties can be exactly substracted out. 

The numerical results presented in the following have been obtained from the direct iteration of the transfer matrix 
equation ft). In Sect. 11111 concerning the end-point transverse fluctuations, we give results for lengths L in the range 
30 < L < 60, with n s (L) independent samples between n s (L = 30) = 35. 10 4 and n s (L = 60) = 23. 10 3 . In Sect. IIVI 
concerning thermodynamic quantities, we have used similar values, as well as smaller sizes to have a bigger range in 
the size L. For instance, the data presented on Fig. Q correspond to lengths L = 6, 12, 24, 48 with respective 
numbers n s (L) of independent samples of order n s (L) = 9 • 10 6 , 9 • 10 6 , 6 • 10 5 , 4 • 10 4 . 

III. STUDY OF SPATIAL PROPERTIES 

The phase transition corresponds to a change in the wandering exponent C characterizing the transverse fluctuations 



< 



i?y (L) > of the end-point of the polymer Q 



< R 2 T {L) > ~ L for T > T c (6) 



< R 2 T {L) > ~ L 2C for T <T C (7) 

L — >oo 

where, as usual, < A > and A denote respectively thermal and disorder averages. The exponent £ ~ 0.59... charac- 
terizes the zero-temperature fixed point 0. On Fig. (JJJ, we show the data for pl{T) = < R^(L) >/L as a function 
of the temperature T for increasing lengths in the range 30 < L < 60. To localize the critical temperature T c , we 
have tried to fit our data according to the following form 



±3M^ =MT)L » +<tl{T)+ <m (8) 

where we have used the scaling relation to — 2( — 1 ~ 0.186. The coefficient ao(T) is found to vanish around 
T c ~ 0.78 — 0.79. A tentative finite-size scaling analysis is shown on Fig. |JSJ|. 
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FIG. 1: (a) p L (T) = < R\{T) >/L as a function of the temperature T for different L = 30 (O) , 36 (□) , 42 (0), 48 (A), 54 
(x), 60 (+) (b) Zoom of the same data in the critical region. 



IV. STUDY OF THERMODYNAMIC PROPERTIES 



For a disordered sample (i) of length L, we note the free energy Fx(L,i), the entropy St(L,i) and the internal 
energy Et{L,i). For each of these observables (X — F,S,E), we will discuss the behavior of the averaged value 
X^ V (L), of the variance AX^ V (L) and of the rescaled probability distribution 



A. Delocalized phase at very high temperature T>T t 




FIG. 2: High temperature phase: (a) Gaussian rescaled probability distribution Gf(x) of the free energy at T = 2 for L — 18, 
30. (b) Non-Gaussian rescaled probability distribution Gs(x) of the entropy at T = 2 for L = 18, 30. 



In the high temperature phase T ^> T c , the free-energy Ft(L, i), the entropy St{L, i) and the energy Et(L, i) of a 
tmple (i) of length L are expected to be 

F T {L,i)=Lf ann + ai (10) 
S T (L,i) =Ls ann + bi (11) 
-&r(£, i) = Le ann + a (12) 
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where dj,&j,Cj are random variables of order O(l). The rescaled distribution J?JJ of the free energy is found to be 
Gaussian (see Fig. 121(a)), whereas the rescaled probability distribution of the entropy is asymmetric (see Fig. 121(b)). 

B. Localized phase at very low temperature T<T C 




FIG. 3: (a) T — fixed point: comparison of the rescaled distributions of the ground state energy Ge {x) and of the free 
energy G F (x) at T = 0.25, 0.4 for L = 12 (b) The rescaled entropy distribution G s (x) for T = 0.25 and L = 18, 30 is found 
to be Gaussian in agreement with the droplet picture. 

Let us first recall some theoretical expectations. At T = 0, the ground state energy of a sample (i) of size L takes 
the scaling form 

E (L,i)=Le +L u Ui (13) 

where m is a random variable of order 0(1). As a consequence, both the correction to extensivity and the fluctuations 
involve the same exponent u> ~ 0.186 

Eo(L,i)-Le =L"ul (14) 

AEo(L) = L" (u> - {ui) 2 ) 1 ' 2 (15) 

Numerical results on the statistical properties of the ground state and low energy excitations are discussed in our 
recent work ■ According to the droplet theory, the whole low temperature phase < T < T c is governed by a zero- 
temperature fixed point. However, many subtleties arise because the temperature is actually 'dangerously irrelevant'. 
The main conclusions of the droplet analysis |]| can be summarized as follows. The scaling 1)1. 'it translates into the 
statistical properties of the free energy, provided one introduces a correlation length £(T) to rescale the length L 

F T (L, i) - Lf ann = + /i(T) (J^j Ut (16) 

where ut is a random variable of order O(l). 

For T low enough (T < 0.4), we find that the width AF(L,T) indeed scales as L 018 . Moreover the rescaled 
distribution © coincides with the rescaled distribution of the ground state energy at T = (see Fig|^a)). 

However, within the droplet theory p|, the random part of order L u of the free-energy (Eq. I16|l is expected to be 
a near cancellation of energy and entropy random contributions of order L 1 / 2 . The argument is that the energy and 
entropy are dominated by small scale contributions of random sign |5j, whereas the free energy is optimized on the 
coarse-grained scale £(T). These predictions for the energy and entropy have been numerically checked in d = 1 and 
d = 2 [a, EH- Moreover, the droplet argument suggests that the entropy can be written as 

St(L, i) - Ls ann = s (T)L + sxi^L^Vi + ... (17) 

where Uj is Gaussian, because the entropy is dominated by a large number of independent local excitations. We show 
on Fig. 01(b) that the rescaled distribution is in good agreement with this argument. 

Then from the thermodynamic relation F = E — TS, the energy is expected to have a more complicated form 
involving terms of order L, L 1 / 2 and L u . 
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C. Critical region 

1. Free- energy scaling at T c 




T T 

FIG. 4: (a) Plot of Y L (T) = F ^^f/ 2 a "" for L = 6 (O), 12 (□), 24 (0), 48 (A) (b) Zoom around the pseudo-crossing point, 
showing a systematic shift with L : the curves for L = 6 and L = 12 cross around T ~ 0.72, whereas the curves for L = 24 and 
L — 48 cross around T ~ 0.755. 

At T c , the wande ring exponent £ is expected to be exactly ( c = 1/2 j24[, and thus the scaling relation ui c = 2( c — 1 
between exponents [23 yields co c = 0. Previous numerical studies and arguments [2TL [2(j have concluded that the 
free-energy is logarithmic at T c with exponent (1/2) 

F Te (L,i) - Lf ann = (hxLf^Wi (18) 

where Wi is a random variable of order 0(1). We show on Fig. 0Ja) that the rescaled variable Yl(T) — Ft ^Lyf™" 
increases with L in the low temperature phase, and decreases with L in the high temperature phase. In the critical 
region, there is a systematic shift in the crossing points of curves when L increases (see the zoom of Fig. Etb)). Since 
our numerical data are not sufficient to analyse quantitatively the L-dependence of this shift, we simply conclude that 
the thermodynamic critical temperature T c is in the region T c > 0.76, and is thus very close to the exact upper bound 
T 2 = 0.79 (Eq. |3. 



2. Rescaled distributions of free energy, energy and entropy in the critical region 

In this section, we compare the rescaled probability distributions of thermodynamic observables in the critical region 
with respect to the low temperature and high temperature corresponding results. 

The results for the free-energy are shown on Figs. [^focusing on the bulk (a) and on the tails (b). The negative tail 
of this probability distribution is of particular interest, since it is linked to the location of the critical temperature 
T c with respect to the upper bound T 2 as explained in details in 0. If this tail is exponential Gf{x) ~ e~ a ^ for 
x — > —00 as on the Cayley tree 0, then T c < T2, whereas if Gf(x) decays more rapidly than exponentially, then 
T c = T<j. For instance, in the mean-field Sherrington-Kirkpatrick model of sp in-g lasses, one has T c — T2 and the 
distribution of the free-energy fluctuations for T > T c is known to be Gaussian |l8l ITgj . Here we obtain numerically 
that Gf{x) is Gaussian in the high temperature phase for T 3> T2 (see Fig. 121(a)). If this Gaussian distribution is 
valid in the whole high temperature phase, then one has the equality T c — T 2 . On the other hand, the strict inequality 
T c < Ti requires that the negative tail becomes exponential in a region above T c . Such a 'transition' in the form of the 
free-energy distribution within a high temperature phase is known to exist for instance in the Random Energy Model 
|20|. Here with our numerical results, as T grows from T2 to T 3> I2, we see a slow crossover from the low-temperature 
tail towards the Gaussian high temperature tail (see Fig. |3J. Since Gf(x) decays in the low-temperature phase as 
e - b \ x \ v w ith 77 = 1/(1 — uj) > 1 according to Zhang argument we see no sign of an exponential decay with 

power 77 = 1 at any temperature. 
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FIG. 5: Free energy distribution : (a) Rescaled probability distribution Gf(x) of the free-energy for L — 12 in the low 
temperature phase T = 0.25 ( thick line)), in the high temperature phase T — 2 (Gaussian, long dashed line) and in the critical 
region T — 0.79 (thin line), (b) Logarithmic plot of the same data: in the critical region, the negative tail is very close to the 
low-temperature tail and far from the Gaussian high temperature tail. 





FIG. 6: (a) Rescaled energy distributions for L — 12 in the low temperature phase T = 0.25 (thick line), in the high temperature 
phase T — 2 (long dashed line) and in the critical region T = 0.79 (thin line), where the distribution is maximally peaked 
(b) Logarithmic plot of the same data : the slower decay in the negative tail is obtained in the critical region and seems to 
correspond to an exponential decay 1hGb(i) ~ — \x\ near T c . 



The rescaled probability distribution for the energy and entropy are shown respectively on Figs. El and [7| focusing 
on the bulk (a) and on the tails (b). Concerning the bulk (see Figs El (a) and (a)), both distributions have a 
non-monotonous behavior as T varies, and become maximally peaked in the critical region 0.7 < T < 0.9. As for the 
negative tails (see Figs El (b) and 13(b)), the slower decay is obtained in the critical region and seems to correspond 
to an exponential decay Ge,s( x ) ~ e _c ' a: '. 
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FIG. 7: (a) Rescaled entropy distributions for L = 18 in the low temperature phase T = 0.4 (thick line), in the high temperature 
phase T — 2 (long dashed line) and in the critical region T = 0.79 (thin line), where the distribution is maximally peaked 
(b) Logarithmic plot of the same data : the slower decay in the negative tail is obtained in the critical region and seems to 
correspond to an exponential decay In Gs (x) ~ — \x\ near T c . 

V. TENTATIVE FINITE SIZE SCALING ANALYSIS ON AVERAGED OBSERVABLES 
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FIG. 8: (a) Plot of the data of Fig [T] for p L (T) = < R 2 L (T) >/L , as a function of the rescaled variable r = (T - T c )L 1/v , 
see Eq. JT^, with T c = 0.79 and ± ~ 0.475. Symbols : L = 30 (O) , 36 (□) , 42 (0), 48 (A), 54 (x), 60 (+) (b) Plot of 
the rescaled energy -SM^ of Eq. (HOJ as a function of the rescaled variable r = (T - T c )L 1/l/ , with T c = 0.79 and ± ~ 0.25. 
Symbols: L = 12 (Q) , 24 (□) , 36 (0), 48 (A), 54 (x), 60 (+). 



Fixing the critical temperature at T c = = 0.79, we have tried to make some finite-size scaling analysis to 
determine the value of the correlation length exponent v, which is related to the specific heat exponent a via the 
hyperscaling relation v = 2 — a. 

We first consider the spatial properties discussed in Sect ion ITTT1 The best rescaling of the data of Fig.l^a) according 
to the finite size scaling form 



< R r (L)> = g( (T _ Te)L i/, 

corresponds to values in the range 0.45 < ~ < 0.5, see Fig. (JSJa). 

We then consider the averaged energy, with the following usual finite-size scaling form 



5E = E^ v {L)-Le ann = L 



(19) 



(20) 
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The best rescaling (see Fig. (0b)) correspond to 0.23 <± < 0.27. 

These two values (as the previous ones found in Ref. 0123) are n °t compatible. As a consequence, we cannot give 
a reasonable quantitative estimate of the critical behavior. In this respect, we recall that an alternative to the usual 
power-law behavior £(T) ~ (T c — T)~ v , has been proposed in ref |l5j ■ namely In £(T) ~ \n 2 (T c — T) + .... Unfortunately, 
our numerical data are not able to test this alternative form without any knowledge of subleading divergent terms. 



VI. CONCLUSION 



In this paper, we have presented numerical results for the directed polymer in a 1 + 3 dimensional random medium. 
Our data concerning the wandering R^(L) of the end-point and the free energy Ft(L) point towards 0.76 < T c < 
T% = 0.79, so our conclusion is that the critical temperature T c is equal or very close to the upper bound Ti derived 
by Derrida and coworkers. We have also presented results for histograms of free-energy, energy and entropy as T 
varies. For T T c , we obtain that the free-energy distribution is Gaussian. For T <C T c , the free-energy distribution 
coincides with the ground state energy distribution, in agreement with the zero-temperature fixed point picture, and 
the entropy fluctuations of order AS ~ L 1 / 2 follow a Gaussian distribution, in agreement with the droplet predictions. 
However, from our various data, we cannot give a reasonable quantitative estimate of the divergence of the correlation 
length £(T) asT^T c 7. 
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